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Аннотация: В статье рассмотрена методика решения изменения мощности внутреннего 
источника тепла с учетом солнечной радиации в многокомпонентных цементных материалах. 
Показана, математическое описание тепловых процессов и предложен метод математического 
моделирования нестационарных полей температуры в многокомпонентных цементных 
материалах. 
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Актуальность. Перспективы существования человечества требуют решения трех 
основных проблем: обеспечение населения пищей, энергией и сохранение природных 
условий, пригодных для жизни. Решение этих задач в значительной степени связано с 
использованием солнечной энергии. 


В условиях дефицита топлива и непрерывного роста цен на него освоение 
колоссальных ресурсов солнечной энергии является одной из важнейших научно- 
технических проблем. 


Исследования, проведенные в последние годы, и опыт по разработке и эксплуатация 
солнечных установок, накопленный в ряде стран, показали возможность и 
экономическую целесообразность более широкого и разнообразного использования 
солнечной энергии уже сейчас, на современном уровень технических возможностей 


[1]. 
Теоретическая часть. 


Мощность объемного источника тепла dy, обусловленного выделением теплоты 
гидратации, изменяется в зависимости от времени, температуры и коэффициента 
излучения в многокомпонентных цементных материалах поли структурного строения 
[2,3,4]. Изменение qy во времени приближенно можно выразить для фиксированной 
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средней температуры кусочно-непрерывной функцией г (рис. 1.а) 


qe >в, (1-1) (1) 


0 


или ломаной 
т 
q~ У (У-Ут-1) ( t-tm)60(t-tm), (2) 
0 
где Oo(T-Tm) - единичная функция Хевисайда, при 
T > Ti OolT-Ti) = 1, прит < п со(т-т) = 0; (3) 
У» - скорость равномерного изменения мощности источника 4 при г-гь, Вт/м?:с; 


Tm - время т-ого изменения мощности источника и скоростей равномерного подъема 
или спада мощности q (рис. 1.6). 


Зная, что: 


1,T>T 
FO) = У qnGo(t-tm), оо) =0 } (4) 


0, T<Tm 


имеем изображение по Лапласу: 


т 


F(S)= У, qmEXP(-stm)/S + Уҹ -K, (5) 


FQ) = У, (У-У) тта): во; 


т К 
Fs= У (Vin-Vin-t)EXP(-stmy/S + Уд: К. (6) 
0 


0 


а) 


Оу, 
BT/M3 
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6) 


ВТАМЗ 


11 Tm Время, г, ч. 


Рис. 1 Аппроксимация действительной кривой изменения мощности внутреннего 
источника q кусочно-непрерывной (а) и ломаной (6) функциями 


Начальное распределение температуры по толщине структурообразующего 
многокомпонентного цементного материала равно to. В начальный момент времени OHO 
помещается в среду с температурой tc > to, которая поддерживается постоянной на 
протяжении всего процесса нагревания Необходимо найти распределение 
температуры по толщине образца и расход тепла в любой момент времени, если 
теплообмен с окружающей средой происходит по закону Ньютона [5,6,7]. Поместим 
начало координат в середине толщины пластины и обозначим через 2/ ее толщину. 
Внутри образца действует источник тепла удельной мощностью Qy, являющийся 
функцией времени [8,9,10]. 


Условие задачи математически может быть сформулировано следующим образом. 
Решить дифференциальное уравнение: 


dt (x, r/d = а drt (x, гах? + 4,(1)/ср; 


(r >0, - 1 <x> 1) (7) 
при условиях 

t (x,0) = t; (8) 

dt (0, т) = 0 (9) 

dx — dt (l, t)/dx + [a// te - fll, t)] = 0. (10) 


В уравнении (7) удельная мощность источника определяется зависимостями (1) или 
(2). 


Применяем к уравнению (7) интегральное преобразование Лапласа [11,12] ... Тогда из 
(5) получим: 
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CP о 


n 1 t 12 exp|— ST,, 
Т, нев bre) o l (11) 


Решение для изображения th(x, 5) при условии (8) 


Т, (x, а а нао (12) 
S 


5 ~ ep 


Постоянную, в дальнейшем A, - найдем из граничного условия (10), которое для 
изображения Т; (х, S) имеет вид 


ЈЕ) (13) 


Удовлетворив решение (11) граничному условию (13), можно определить постоянную 
«А» [13,14]... Из зависимости (12) получаем 


5 1 
T, (e) онь ; (14) 
a \а 


а прих =/ 


1 
T, (1,5) = аа Е (15) 
а Va 


Подставляя значение производной Т; '(/, 5) в условие (13), получим 


5 S а 5 A\t.-t 1 < exp(—sz,, ) 
A h,|—l-—ch,|-1| = c__0 _ ; (16 
a j T | al 5 252-4 5? | oo 


Й a Ta) 


t, — 1, ср o 


(17) 
Jaf 1+ а [fè i) а сто) 
Тогда решение (12) примет вид: 
5 
t —t \ch —х 
t 1 expl- sz, (1, 0 E 
T,(x,s)- я = 252.4, E Я 2 
: б T] 
а AN\a a 
(18) 


1“ exp|— ST n 
У 4, | 2 ) 


ср 5 
fit Ef 
a AN\a a 


или, переходя к температуре t (x, T), получим 
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m 2 
g 907) Brey r, )P |152 |- 


Е 2 Om P B. 
t. — to 0 i (19) 
a р: А, . ‘совр, техр(- wer) 
n=l ш! 


= Ae- Е =. 


Если рассматривается не мгновенный подъем температуры на границе изделия, а 
постепенный, то для первого периода подъема температуры со скоростью 


V= (tinax-to) /To 


имеем: 


оў ў ye A . ехр(- i) + 


n=1 n 


t(x,7) = Е 
(21) 


т 


И ЕБЕ ИИ 


n=l п 
а динамика температуры в точке х = 0 (посередине образца) 
(22) 


Для второго периода (изотермическая выдержка): 


1 1 
t (0, т) = Үт. È 2 [(1-Q1, г)-( 1-Ф1, 0)] + 2 У (ke N-1) (= Tk) [1-ФКт- T%)] (23) 
0 


Обсуждение: Решение получены для двух периодов гелиотеплохимической обработки: 
равномерного плавного подъема температуры среды в гелиотехнологической камере и 
поддержание ее на определенном уровне. 


Вышеуказанные аналитические решения позволяют качественно оценить учет 
коэффициента луче поглощения солнечной радиации в развитие поля температур в 
изделии с выделением теплоты при протекании экзотермической реакции твердения 
вяжущего. 


Выводы: на основе принципов системного анализа разработано математическое 
описание тепловых процессов и предложен метод математического моделирования 
нестационарных полей температуры в многокомпонентные цементные материалы. 
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